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Abstract 

-xj- : 

. Hamiltonian structures for non-linear hydrodynamics of incommensu- 

1 rate intergrowth compounds (IIC) and quasicrystals (IQ) are constructed. 

We discuss also the way to account for internal friction of phason nature. 
We show that the existence of a self-force in IIC and IQ is not only matter 
■ of constitutive issues, rather it is related with questions of SO (3) invari- 

' ance. The covariant mechanics of discontinuity surfaces in quasiperiodic 

structures is also analyzed. The attention is mainly focused on the inter- 
action between 'diffuse' grain boundaries and sharp discontinuity (moving 
possibly) surfaces, 
•i— I . 

' To T. Y. for her moral vigour to oppose racialism of stupid people and 

, to face difficulties. 

1 Introduction 

Quasiperiodic metallic alloys display two types of low energy excitations in the 
hydrodynamic range: the standard phonon modes associated with the congruent 
distortions between neighboring material elements, and phason modes (Gold- 
stone degrees of freedom) due to local rearrangements of atomic clusters. The 
former modes are represented by the standard displacement field u, while the 
latter ones by the phason vector field w. 

Here our attention is focused on incommensurate intergrowth compounds 
(IIC) and proper icosahedral quasicrystals (IQ). For IIC, the phason displace- 
ment field describes collective atomic modes associated with relative displace- 
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ments between incommensurate sublattices determining quasipcriodicity. Basi- 
cally, IIC can be considered as the result of modifications of periodic structures. 
On the contrary as shown by the 1984 experiments of D. Shechtman, I. Bleck, 
D. Gratias and J. W. Cahn [S], there exist intrinsically quasiperiodic crystals 
(IQ) that are not modulation of multiply twinned periodic crystalline struc- 
tures or the composition of more species of them. They admit icosahedral phases 
with consequent long-range orientational order and absence of translational one. 
They violate the crystallographic restriction prescribing that a crystal displays 
periodicity and "cannot have certain forbidden symmetries, such as fivefold ro- 
tation" [L] but do not constitute a new state of matter. For IQ, the phason 
displacements describes substructural changes of diffusive nature [RoLo]: (i) 
collective atomic modes and (ii) tunneling of atoms below energetic barriers 
separating places at a distance lesser than the atomic diameter. 

Phason activity thus exists in IIC and IQ; however, the energetic landscape 
is different. In the diffraction scenarios obtained by x-ray scattering experi- 
ments, diffuse scattering is registered around Bragg peaks in both cases. In 
principle, such a scattering could be represented through the singular part of 
certain autocorrelation measures [BH]. However, for IIC there are six sound- like 
branches while in the case of IQ just three sound-like branches appear. As a 
consequence, kinetic energy can be attributed to the phason activity in the case 
of IIC, while for IQ there may be a sort of internal friction leading to viscous-like 
evolution because phason activity displays diffusive nature. Since in the case of 
IIC, at each point X, the vector w = w (X) represents the relative displacement 
of incommensurate sublattices, it is a measure of deformation and does enter 
the structure of the free energy together with its gradient Vw. It docs not 
happen to IQ where just Vw appears in the list of constitutive entries of the 
energy as a representative of phason behavior. The explicit dependence on Vw 
varies according to the circumstance that IQ is in a 'locked phase' (i.e. a phase 
without phason contribution of Debye- Waller type) or in an 'unlocked phase' 
(see [JS]). 

Taking into account analogies and differences, we present here a Hamiltonian 
formalism for non- linear hydrodynamics of IIC and IQ. We follow strictly the 
general framework of multifield theories [GC], [C], [M], [CM] unifying a wide 
class of models of condensed matter with complex substructural morphology. 
In discussing the non-linear elasticity of quasiperiodic structures, we follow in 
particular the general results in [CM] dealing with field theories in which the 
Hamiltonian accounts for order parameters taking values on an abstract man- 
ifold. Here, at each point X, our order parameter is the phason displacement 
w = w (X) , belonging to some copy of the translation space of the three- 
dimensional Euclidean point space £ 3 . 

The mechanics of quasiperiodic crystalline structures in the hydrodynamic 
range, the one of elasticity and plasticity, has been discussed variously in scien- 
tific literature. The attention has been focused primarily on quasicrystals rather 
than on IIC starting from the work of T. C. Lubesky, S. Ramaswamy, J. Toner 
[LRT] on (see [HWD], [RT], [DP], [JS], [RL], [RoLo], [DYHW], [MSA]). 

Here, we analyze the matter. Accomplishments are briefly summarized be- 
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low. 



(i) Hamiltonian structures for both IIC and IQ are constructed. The covariance 

of the balance of phason interactions follow. A balance of interactions 
occurring when defects in the quasiperiodic structure are permuted is then 
deduced. 

(ii) It is possible to show that the existence of a self-force of phason nature 

within each material element is a consequence of SO (3) invariance. It 
is characterized by constitutive instances: in the case of IQ it vanishes 
at thermodynamical equilibrium. A non-standard integral balance of mo- 
ments follows naturally from SO (3) invariance. It is different from the one 
postulated commonly (see [HWD]) and gives rise (by localization) to the 
pointwise balance of phason interactions that is not necessarily associated 
with an integral balance of phason momentum. 

(iii) We discuss also non conservative issues to account for the internal friction 
in IQ. Our treatment contains the 'minimal model' presented by S. B. 
Rochal and V. L. Lorman in [RoLo]. 

(iv) Finally we describe the influence of phason activity on the evolution of 
discontinuity surfaces in quasiperiodic crystalline structures. We show 
covariance of the balances of phonon and phason interactions at the dis- 
continuity surfaces. 

For the sake of simplicity our treatment deals only with isothermal processes. 

Algorithms for analyzing numerically specific cases follow naturally. 

Unaspectcd phenomena may be evidenced. For example, in the linear case, 
due to uncertainties in the experimental evaluation of phonon-phason coupling 
coefficient in icosahedral quasicrystals, the combined use of Monte-Carlo and 
finite clement techniques allows us to put in evidence the possible stochastic 
clustering of phonon and phason modes around macroscopic defects like cracks. 
The relevant results will be presented in a forthcoming paper with M. Gioffre 
and F. L. Stazi. 

Some notations. For any pair of vector spaces A and B (with duals A* and 
B*), Horn (A, B) is the space of linear maps from A to B. For any manifold 
M, T m M is the tangent space of M at m e M, while T^M the relevant cotan- 
gent space. Moreover, Aut (A) indicates the space of automorphisms of A. We 
will make use of two different regular bounded regions of the three-dimensional 
Euclidean point space £ 3 , namely Bo and B, and of two different copies V u and 
Y w of the translation space of £ 3 (wc may also identify them as copies of R 3 ). 
Capital letters A, B, C... used as indices denote coordinates in Bo, while k... 
coordinates in B. The differential operators Div and V indicate respectively 
divergence and gradient calculated with respect to coordinates in Bo while div 
and grad are their counterparts with respect to coordinates in B. The super- 
script T means transposition. The symbol d y means partial derivative with 
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respect to the entry "y" . We indicate with the term part any subset of Bq with 
non- vanishing volume and the same regularity properties of Bo- Let £ be any 
smooth surface in Bo oriented by the normal m at each point, for any field e (•) 
defined on Bq and diffcrcntiable there, we indicate with Vs its surface gradient 
along £, namely Vse (X) = Ve (X) (I — m® m), with I the second-order unit 
tensor. The trace of V^e is the surface divergence of e, namely DivY,e. Other 
notations will be explained later. 

2 Configurations, observers and relabeling 

2.1 Configurations 

Our analysis deals with a quasiperiodic crystalline body occupying in its refer- 
ence place a regular 1 region Bq of the three-dimensional Euclidean point space 
£ 3 . A generic point X £ B is identified with the centre of mass of a crys- 
talline cell (which is the characteristic material element) that one may imagine 
collapsed at X in a coarse grained representation of the quasiperiodic structure. 

A standard deformation of the body is represented by a sufficiently smooth 
injective mapping Bo 9 X x = x(X) £ £ 3 . The current place B = x (Bo) 
of the body is a regular region too. The placement map x is also orientation 
preserving: at each X its gradient F = Vx, i.e. the value of the field Bo 3 

X^F = F(X)e Horn (T^B , T X B), has positive determinant. 

Let g be the spatial metric in B and 7 the metric in Bo- The linear operator 
F T F = C £ Hom{Tx_B ,T*B) is the pull-back at X of g through x, i.e., in 
coordinates, Cab = F^dij^B- Then, the difference (C — 7) is twice the non- 
linear deformation tensor E. 

If we consider each material element as a perfect crystalline cell, during a 
motion Bo x [0, t\ 9 (X,i) x = x (X,t) £ £ 3 , the standard displacement field 
u = u (X,t) = x (X,i) — X £ Y w is the descriptor of phonon degrees of freedom. 
When the material element undergoes at least one of the substructural changes 
(i) and (ii) described above, namely collective atomic modes or tunneling of 
atoms, a sort of internal shift occurs and is represented by a vector w so that 
we have a vector field Bo 9 X w = w (X) £ Y w that we presume sufficiently 

smooth over the body. During a motion, we then have Bo x [0, t\ 9 (X,i) ^—> 
w = w (X,t) £ V w , with a slight abuse of notation. 

From the point of view of the general setting of multifield theories, the copy 
V„, of the translation space V over £ 3 , containing w, plays the role of the mani- 
fold of morphological descriptors (order parameters) of the material substructure 

n [m]. 

One may consider (see [HWD]) a global displacement u belonging to 'Vu'S^w 
By indicating with x' = x' (X) the point given byx' = X + u = X + u + w, and 
with F' the gradient Vx' (X), we get additive and multiplicative decompositions 
given respectively by F' = F + Vw and F' = F ph F, with F ph = I + (Vw) F" 1 . 

1 Bo is regular in the sense of D-regions defined in [D] . 
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Really, since x is one-to-one, one may construct a representation of w on the 
'apparent' current place B of the body. By indicating with w a = w a (x) 
the image of w attached at x = x (X) e B, we get w a = w o x -1 so that 
(Vw)F _1 = gradw a . As a consequence, the interpretation of the multiplica- 
tive decomposition of F may be the following: we may deform first the body 
at a coarse grained level maintaining frozen phason activity, then we may al- 
low collective atomic modes to develop. In other words, by indicating with 5c ph 
the mapping x ph = x' o xT 1 , we see that F ph = I + gradw a is the gradient 
of deformation from B to Z ph (B), namely there is a piecewise continuous map 

F ph such that B 3 X £^ F ph = F ph (X) e Horn (T x B,T x ,Z ph (B)) . The 
map x' describes the circumstance that collective atomic modes or tunneling of 
atoms occurring within each crystalline cell may shift the centre of mass of the 
crystalline cell itself from its current (in certain sense 'apparent') place x. 

Finally, we indicate with x = ^x (X,t) and w = (X,i) rates in the 
reference description and might use also u = ^u(X,i) instead of x to put in 
evidence the role of phonon and phason degrees of freedom. 

If we restrict our attention to infinitesimal deformation regime in which 
B can be 'confused' with Bo in the sense that x « u at each X, in addition 
to the standard compatibility condition curl curl symVu = 0, we get also a 
phason compatibility condition curl curl Vw = that would imply eventually 
an energetic contribution of phason spin. 

2.2 Observers and relabeling 

For the mechanics of quasicrystals the definition of the concept of observer 
follows general issues of the mechanics of complex materials (see [M04]) involving 
in such a definition the representation of all geometrical environments necessary 
to the description of the material morphology. 

Three sets enter in fact the geometrical picture of a quasi-periodic crystalline 
body: the point space £ 3 (i.e. the standard ambient space), the translation space 
V„, (containing phason degrees of freedom) and the interval of time [0,t\. An 
observer O is then a representation of £ 3 , Y w and [0,t\. 

We consider also relabeling of material elements in Bo, simulating a redistri- 
bution of possible defects. 

Relabeling. Formally, a 'permutation of inhomogeneities' in Bo is described 
by the action of the special group of isocoric diffeomorphisms SDiff on Bo- So 
that we have a map 

• R+ 3 si i — > f s \ £ SDiff (B ), with fo 1 the identity. 

At each s\ we get X i — ► f Si (X), with Divf}^ (X) = 0, where the prime 
denotes differentiation with respect to the parameter s\. We put fg' (X) = ro. 

Changes of observers. We consider observers agreeing about the measure 
of time so that a generic change of observer involves just a couple of transfor- 
mations: one of the ambient space £ 3 , the other of V w . They are described by 
the parametrized families of mappings defined below. 
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• K+ 3 s 2 ' — > f 2 2 G {£ 3 ) , with f 2 the identity. We put f 2 ' (X) = v. 

• A Lie group G, with Lie algebra g, acts over V w . If £ e g, its action 
over w EW W is indicated with £ Vm (w). By indicating with w s the value 
of w after the action 2 of g e G, if we consider a one-parameter smooth 
curve M+ 3 S3 1 — ► g S3 € G over G such that £ = | S3= o and its 
corresponding orbit Si 1 — ► w Ss3 over V w , starting from a given w, we 
have£ Vw (w) = ^w 9s3 | S3=0 . 

3 Lagrangian structures for phonon-phason elas- 
ticity 

Up to this point just geometry has been involved. In constructing a mechani- 
cal model of a body after the description of its morphology, one discusses the 
representation of interactions and their balance first, then the explicit represen- 
tation of constitutive relations. The two issues are essentially separated. The 
representation of interactions by means of appropriate vectors or higher order 
tensors is a consequence of the essential geometrical description of the body (in- 
teractions are in fact entities power conjugated with the rates of morphological 
descriptors) and the balance is independent of the constitutive nature of the 
material. 

When we develop Lagrangian and Hamiltonian formalisms as below, in in- 
troducing the Lagrangian density just after geometrical issues, we put on the 
same ground the representation of interactions and constitutive issues because 
they are mixed in the variational description. 

Let us consider a fiber bundle 

n:y^B x [0,i] (1) 

such that 7T" 1 (X,t) = £ 3 xY w is the prototype fiber. A generic section r\ e T (y) 
is then a mapping 77 : Bo x [0, t\ — ► y such that 77 (X,t) = (X, t, x, w) with x 
and w in the fiber n^ 1 (X,t). If sufficient smoothness for sections is allowed, 
the first jet bundle J x y over y is such that 

J'y 3 j 1 fa) (X,i) = (X, t, x, x, F, w, w, Vw) . (2) 

In the conservative case we presume that the canonical Lagrangian 3 + 1 
form 

L : J x y — > A 3+1 (Bo x [0, t\) (3) 
admits a sufficiently smooth density C such that 

L (j 1 (77) (X,i)) = C (X, t, x, x, F, w, w, Vw) d (vol) Adt. (4) 

2 It is not essential to render precise if the action is from the left or from the right. 
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with C defined by 

1 2 1- 2 
£(X ) i,x,x,F,w,w,Vw) = -p |x| +-p|w| - 

-p e(X,F,w,Vw) -p w(x), (5) 

where p is the referential mass density (conserved during the motion), p an 
inertia coefficient for possible phason kinetics (see [HDW], [RL]), e the elastic 
energy density and w the density of the potential of external actions, all per unit 
mass. Here we do not consider possible bulk external direct actions on phason 
changes. 

We then evaluate the variation of the total Lagrangian L (Bo) given by 
L (Bq) — /g oX [Qj] Cd(vol) A dt and we may find at least one section (with the 
properties of x and w) satisfying Euler-Lagrange equations for L (Bo), namely 

d^C = d x C- DivdpC, (6) 



<9w£ = <9 W £ — Divdy^C. (7) 

Definition 1 (invariance of C). C is invariant with respect to the action of 
fi, ft 2 and G if 

C (X, x, x, F, w, w, Vw) = 

= C (f\f 2 , (gradi 2 ) x, (gradi 2 ) F (Vf 1 )" 1 , w ff , w s , (Vw,) (Vf 1 )^) . (8) 

where we indicate with f 1 , f 2 and w g the values (X), f 2 2 (x), w 9s (X). 

Let Q and 5" be scalar and vector densities given respectively by 

Q = diC- (v - Fro) + d*£- (£ Vra (w) - (Vw) to) , (9) 

d = Cxv+ (d F C) T (v - Fro) + (d Vvi C) T (£ Vw (w) - (Vw) n) . (10) 

Theorem 1. If the Lagrangian density C is invariant under fi, f 2 2 and G, 
then 

Q + Div$ = 0. (11) 

Theorem above is a version for quasiperiodic bodies of Noether theorem. A 
generalization of it for multifield theories that involve order parameters belong- 
ing to abstract manifolds is proven in [CM]. 

Corollary 1. If f 2 2 alone acts on C leaving v arbitrary, from 1|11[) we get in 

covariant way the balance of phonon interaction (standard Cauchy's balance of 
momentum) 

p x =p Q h+DivP, (12) 



7 



where P = —&fC is the first Piola-Kirchhoff stress and b = d x £ the vector of 
body forces. 

At each X in Bq, P maps linearly normals to surfaces through X into ten- 
sions at x in B. There is then a map P such that Bq G X i — ► P = P (X) E 
Hom(T^B ,T*B). 

Corollary 2. If G arbitrary acts alone on C, from we obtain the balance 
of phason interactions 

pw = — z + DivS, (13) 

in covariant way where S = — dy w C represents phason stress due to the rel- 
ative influence of the phason activity between neighboring material elements; 
z = —p d w e (self-force) describes self-interactions of phason nature within each 
material element. 

At each X, S e Hom(T£B ,T£V w ) and z e T^¥ w . The phason stress 
S maps linearly normals to surfaces through X in Bq into tensions of phason 
nature, i.e. elements of T^,V W ~ K 3 . 

Corollary 3. Let G = SO (3) and, for any element qA of its Lie algebra, 
be such that v = qA (x — xo) with xo a fixed point in space. If C is independent 
of x and only the special choices of ff and G just defined act on £, one gets 
from ljTT|) 

skw (<9 F eF T + w <£> <9 w e + (Vw) <9 Vw e T ) = 0, (14) 
where skw (•) extracts the skew-symmetric part of its argument. 

Corollary 4. If alone acts on C, with ro arbitrary, from (|llfl one gets 
(F T diC + Vv T do£) - Div MP-i (p |x| 2 + p |w| 2 ) - dx£ = (15) 

where P = p el — F T P — Vw T 5, with I the second order unit tensor, is a 
generalized Eshelby tensor accounting for phason activity (a special case of the 
general one obtained in [M]). 

Corollary 5. Let G — SO (3) and, for any element qA of its Lie algebra, 
fg is such that w =q A (X — Xo) with Xo a fixed point in Bq. If the body is 
homogeneous, and only the special choices of f 2 2 and G just defined act on C, 
¥ is symmetric. 

Remark 1 (reduction to IQ). In the case of IQ, the elastic potential e does 
not depend on w and inertial effects associated with phason modes are absent. 
So, in the purely conservative case the balance of phason interactions and l|14fl 
become 

DivS = 0, skw (<9 F eF T + (Vw) <9 Vw e T ) = 0, (16) 
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respectively. 



Remark 2. Notice that for fixed w and Vw a standard result of non-linear 
elasticity must hold: namely e (■, w,Vw) cannot be convex in F for reasons of 
SO (3) invariance when large deformations occur. The same property needs to 
be satisfied by e (•, Vw) in the case of IQ. 

3.1 Universal phonon-phason changes in quasicrystals 

For an icosahedral quasicrystal (IQ) we say that the deformation is affine when 
F' does not depend on X. 

Moreover, if we consider deformations that can be controllable just by ap- 
plied macroscopic tractions, excluding in this way body forces, we call universal 
all deformations that can occur in these conditions for all bodies in a given class. 

Theorem 2. Let the mappings (F, Vw) i — > P = P (F,Vw) and (F,Vw) i — ► 

S = S (F, Vw) admit bounded partial derivatives with respect to their entries 
and, at each X, 

det (?i %?s)*«- '"> 

All universal static deformations of homogeneous (purely) elastic quasicrystals 
satisfying the restriction l|17|) are affine. 

Such a theorem is in a certain sense a middle generalization of a standard 
result in non-linear elasticity theory of simple bodies (see [A], p. 506). Here 
the difference relies upon the circumstance that phason degrees of freedom are 
involved. 

To prove it, first recall that in absence of body forces and in conditions of 
homogeneity, for IQ the equilibrium equations read 

DivP = 0, DivS = 0. (18) 

We have also 

DivP = (d F P) VF + (<9 Vw P) VVw = 0, (19) 

DivS = (9 W 5) VF + (<9 V wS) VVw = 0, (20) 

i.e. 

AiVF + A 2 VVw = 0, (21) 

A 3 VF + A 4 VVw = 0, (22) 

with A^'s fourth-order tensors that are arbitrary because the explicit form of 
the mappings P and S is not specified. They are also constant because the 
material is homogeneous. Consequently, thanks to (|T7|) the solution to l(2*T|) and 
l|22|) provides VF = and VVw = 0, that is F and Vw must be constant. As a 
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consequence, thanks to the additive decomposition F' = F + Vw, F' is constant 
as well. 

3.2 Mutations of material metric 

The extended Eshelby tensor P = p el — F T P — Vw T £, accounting for phason 
effects, enters the picture of the interactions involved in mutations (such as 
evolution of defects, interfaces etc.) [M] . These mutations may be represented by 
means of 'alterations' of the geometrical structure of Bo which, on the contrary, 
would remain fixed once and for all. In particular, we focus here our attention 
just on mutations that may involve changes in the material metric 7 defined 
on Bo and energy associated with them. Such a kind of situation may occur 
in plastic flows (see [CM98]) or in the alteration of possible pre-stressed states 
(see, e.g., [Ml]). From now on we assume in this section p = 1 for notational 
convenience. 

We then consider (with some slight abuse of notation) a density of elastic 
energy of the form 

e = e(7,F,w,Vw) (23) 

in which we express explicitly the presence of 7 and require that e is invariant 
under the action of the group of point-valued diffeomorphisms defined on B and 
altering it. 

In other words, we require that e is invariant under virtual superposition 
of sufficiently smooth deformations altering the reference configuration. Notice 
that this kind of request of invariance is more than the request of invariance 
with respect to relabeling because here it is not required that diffeomorphisms 
involved are isocoric. 

To this end, we then consider a one-parameter family h s of sufficiently 
smooth point-valued diffeomorphisms defined over Bo and indicate with u and 
H s the derivatives | s= o and Vh s , respectively. After the action of h s , the 
density e changes as 

e^e hs =(dctH s )p:(H s T T 7 H; 1 ,FH s - 1 ,w,(Vw)H; 1 ). (24) 

Theorem 3. If for IIC and IQ the energy density e depends on the metric 
7 in Bo and is invariant in the sense defined above, one gets 

7?i = 2{d^e) AC lcB . (25) 

The proof relies upon the circumstance that the requirement of invariance 
of e under the action of h s implies j^eh s \ s =o ■ It coincides with 

e (7, F, w,Vw) trVu - d y e ■ ((Vu) T 7 + 7 (Vu)) - 

-d F e ■ FVu - 9 Vw e • (Vw) Vu = (26) 
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since gjH s 1 | s=0 = — Vu. As a consequence of the symmetry of 7, we then get 
(el - F T P - Vw T S - 2 {dye) 7) • Vu = 0, (27) 

which implies Q2fijl. thanks to the arbitrariness of Vu. 
Remark 3. Notice that 

mAD mA BD (Pi \ AC BD r, fa \AD , oq n. 

P = P s7 = 2 {dye) 7 CB 7 = 2 (<9 7 e) (28) 
is symmetric since 7 does. 

4 Elementary Hamiltonian structures 

Hamiltonian structures follow in a natural way from the Lagrangian representa- 
tion described so far. Let in fact p and \x be respectively the canonical momen- 
tum and the canonical phason momentum defined by p —d^C and fi —d^,C. 
The Hamiltonian density TL is then given by 

H (X, x, p, F, w, n,Vv>) = p x + n ■ w - C (X, x, x, F, w, w,Vw) . (29) 

In terms of partial derivatives of TC, the balances I© and J7J) can be written 

as 

p = —d x H + Divd-pH.^ 

x = d p H, (30) 
(1 = -d^H + Divdx/ W H, 

W = dfjTi., (31) 

which are Hamilton equations for IIC. In the case of IQ, w disappears in the 
list of entries of H and l|31|) reduces to 

fi = Divd^wTC, 

w = df/H. (32) 
Both in the case of IIC and IQ, general boundary conditions of the type 

x(X)=x ond {x) B , (33) 

dpHn = t on d {t) B , (34) 

w(X)=w ond {vr) B Q , (35) 
d Vw Hn=t ond (t) Bo; (36) 
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hold, where x, t, w and t are prescribed on the relevant parts d^Bo of the 
boundary, chosen to be such that d^B nd^B = with CI (dB ) = CI (d^B U 0(* 
and d^B n d^B = with CI {dB Q ) = CI (d^B Q U d^B ), where CI indi- 
cates closure and n is the outward unit normal to dBo at all points in which it 
is well defined. 

Hamilton equations above are special cases of the ones discussed in [CM], 
where general order parameter fields taking values on an abstract manifold are 
accounted for. 

It is rather difficult to imagine a loading device prescribing phason tractions 
t at the boundary. Problems with traction data might thus involve the existence 
of at least one surface density U (w) such that t = p^d^U if not another density 
of the type U (x) with t = p Q d x U. In this way one considers the external 
boundary as a structured surface enveloping the body. 

In this case, the Hamiltonian H of the whole body is then given by 

H (x, p, w, (jt) — J H (X, x, p, w, fi) d (vol) — 
Jb„ 

(U(x)-U(w))d(area). (37) 

/<9( 2 >Z3o 

where O^Bq = d^Bolid^Bo- Notice that we write H (X, x, p, w, /i) instead of 
H (X, x, p, F, w, /i,Vw) because below we consider directly variational deriva- 
tives. 

Theorem 4. The canonical Hamilton equation 

F = {F,H} (38) 

is equivalent to the Hamiltonian system of balance equations <|30|) - (|31[1 for a 
quasiperiodic body where F is any functional of the type J B f (X, x, p, w, fi), 
with f a sufficiently smooth scalar density, and the bracket {•, •} is given by 



J Bo \5x Sp Sx Sp 

f f^L^H 

J aw Bo \ Syi S P 



d (vol) + 
8U Sf 



SMflo ~ ■ ^ \dWBo )d(area) 



sf_ sn l sn sf 

d a) Ba \Sw Sfj, 



sn sf , \ 

aWB„ - fa ■ 1- \dWB J d ( area ) 



sf sn sn sf\ 7/ 

where the variational derivative ^ is obtained fixing p and allowing x to vary; 
an analogous meaning is valid for the variational derivative with respect to the 
phason degree of freedom. 

The proof follows by direct calculation (see a more general version of this 
theorem in [CM]). 
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Remark 4. The bracket {•, •} is bilinear, skew symmetric and satisfies 
Jacobi identity. 

Remark 5. For F = H, 1|38|) coincides with the equation of conservation of 
energy. 

Let us consider a boundary value problem in which traction data are not 
prescribed at the boundary where just conditions like l|33|) and (|35|l hold. In 
this case reduces to 

{F,H}= f {f,H} P d(vol), (40) 

J Bo 

where 

ttv\ ( 5 -l 5 Ji_ 5 2i 5 J\<(^_ ^H_sn sn 

lJ ' ip ~ \Sx ' Sp 6x ' Sp) + \Sw ' Sw'SfiJ- [ ' 

It is matter of simple calculation to verify that {-,-}p is bilinear, skew- 
symmetric, satisfies Jacobi identity and also Leibniz identity. Then, {•, - } p 
induces relevant Poisson structures. 

A purely spatial representation of all structures above described is available 
with all fields defined over B rather than Bq. Then Bq does not appear more 
neither as a reference place nor as a model of paragon for quantities involved in 
the mechanical model. In a spatial representation we start assuming a structure 
for Hamiltonian density of the form 

^(x,p,g,w a ,/x a , ff radw Q ) , (42) 

where now, at a given t, (J, a is the canonical momentum conjugated with w a and 
TL depends also on the spatial metric g rather than the gradient of deformation 
F because no reference to Bq is made. In this case, Hamilton equations read for 
IIC 

p = -d±H + div (^2d s H - (gradw) T <9 ffra dw?^ , 

x = d p H, (43) 

A = -d^K + divdgradw'H, 

w = d„H, (44) 
with the obvious reduction to IQ. Here, the Cauchy stress a is then given by 
er = 2<9 g H - (gradw) T d gradvr H. (45) 

The term (gradw) T 9 ffra dwW rules an exchange of energy between the gross 
scale of macroscopic deformation and the finer scale of phason changes, and vice 
versa. An analogous phenomenon occurs in complex fluids where topological 
transitions in the flows may be generated by this type of energy exchanges 
[M03]. 



13 



5 Surfaces of discontinuity: conservative behav- 
ior 



Surfaces across which some quantities undergo finite jumps may occur in IIC 
and IQ. They may be shock or acceleration waves, dislocations, closed cracks 
and so on (see also remarks in [ML]). To describe their behavior, say rules for 
their potential evolution, one needs not only to account for their geometry and 
the action of standard interactions due to deformation processes, but also to 
phason interactions. 

Below we consider discontinuity surfaces endowed with own surface energy 
which models the physical circumstance that interfaces arc often regions in a 
mctastable state with an high concentration of energy [ML] . We allow not only 
discontinuities of the standard gradient of deformation F across the surface but 
also jumps of w and its gradient. Really, one may argue that the presence of 
Vw in the list of entries of the energy takes into account in a regularized way 
the possible presence of grain boundaries. This is true when the grain boundary 
is between two regions with constant phason activity. However, in presence of 
defects or in presence of subgrains containing 'worms' (i.e. the topological al- 
terations of lattices assuring quasiperiodicity), we may have interaction between 
diffuse interfaces and sharp discontinuity surfaces. This is exactly the situation 
that we are analyzing here. 



5.1 A single discontinuity surface in £> 

Let E be a single surface coinciding with {X G clBo, g (X) = 0}, where g is 
a smooth function (smoothness chosen for the sake of simplicity) with non- 
singular gradient. It is oriented by the normal vector field m = m (X) = 
V 5 (X)/|V.g(X)|. 

Let X i — ► a = a (X) be a field taking values in a linear space and suffering 
bounded discontinuities across S. For e > we indicate with a ± the limits 
lim £ ^o a(X±em). The jump [a] of a across £ is defined by [a] = a+ — a~ 

while the average (a) by 2 (a) = a + + a - . If fields a\ and a2 have the same 
properties of a we have [0102] = [ai] (02) + (ai) [02] with the product a±a2 
defined in some way assuring distributivity. 

£ is coherent when at each X one gets [F] (I — m ® m) = 0. 

For any 'virtual' motion to £ prescribed by means of a vector field 

£ 3 xA ro = w (X) e M 3 (46) 

with normal component U = vj ■ m and assume that the velocity x may suffer 
bounded jumps across £, we get the condition [x] = — U [F] m. 

At each X G X we define the surface deformation gradient F as (F) (I — m <g) m) £ 
Horn (TxS,T x S) and indicate with N the projection over S of the average of 
Vw, namely (Vw) (I - m ® m) 6 Horn (T X S, T W Y W ). 
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5.2 Phonon and phason surface measures of interaction 
and their balance 

We consider E endowed with a surface energy density <j) assumed to be suffi- 
ciently smooth and given by 

(m,F, (w) ,N) ^ 4> = 4> (m,F, (w) ,N) . (47) 

The dependence of on m accounts for anisotropy of E. 
We require the invariance of </> with respect to 

(i) changes of observers and 

(ii) relabeling of E. 

As discussed above, changes of observers are characterized by the action of 
the group of automorphisms of £ 3 and of a generic Lie group over V„, . However, 
the definition of fj- needs to be modified in order to describe the relabeling of 
E in addition to the overall relabeling of £> . 

We should then consider time-parametrized families s\ i — > f} of elements 
of SDiff (Bo) characterized by the properties listed below (see [dFM]). 

1. The map si i — ► f* satisfies Al. Moreover, the field Bo X i — ► n> = 
It) (X) = fg (X) is at least of class C 1 (Bo), then across and along E. 

2. Each fj- preserves the elements of area of E. Namely, if dA is the element 
of area of E in 8 0l dA = f** o dA, where the asterisk indicates push 
forward. 

3. (Vtu)m = 0. 

4. Vsv m = 0, with v TO = tn ■ m. 

Definition 2 (invariance of 4>). A surface energy density <f> is invariant with 
respect to the action of , and G if 

~4> (m,F, (w) ,N) = ^Vf 1T m, (gradrf 2 ) ¥ (vf 1 ) , (w) g ,N 9 (vf : ) , 

(48) 

for any g £ G and si, s 2 G M + , where N s = (Vw s ) (I — m ® m) and we have 
used notations common to Definition 1. 

Let X be a sufficiently smooth vector density defined over E by 

X = -0riw+(cW) T (v - (F) w) + (ch<j)) T (£ v ((w» - (Vw) to) -(d m <j) ® m) to. 

(49) 
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Theorem 5. Let S be a structured surface with surface energy (j>. Let us 
assume 



-j- / Qd (vol) + f $ ■ nd (area) + f X ■ nd (length) = 



(50) 



for any part bs of Bq crossing E. If C and <f> are invariant with respect to f * , 
fg 2 and G, covariant pointwise balances across £ follow as in the list below. 

1. TTie action of f% alone implies the interfacial balance of standard inter- 



actions 



where T = — 



[P] m + £>w E T = -p [x] U, (51) 
G Horn (TxS, T^B) is i/ie surface Piola-Kirchhoff stress. 



2. TTie action of G alone implies the interfacial balance of substructural in- 
teractions 

[S] m + Div^S -$ = -p [w] [/, (52) 

w/iere § = — <9n0 £ flom (^xS, T^V W ) is the surface microstress and 3 = 
d w (f) £ T^¥ w the surface self- force. 

3. The action of fj alone implies the interfacial configurational balance along 
the normal m in absence of dissipative forces driving S, namely 

m-[P]m + C tan -L + £)w E c = 



pU 



(Vw) T w 



m+-p 



w 



where 



|Fm| 



C tan = </>n - F T T - N T S 
is a generalized version of the surface Eshelby stress 

c = -d m <p - T T (F) m - § T (Vw) m 

is a surface shear. 



(53) 
(54) 

(55) 



An analogous theorem can be found in [dFM]. However, though there order 
parameters taking values in an abstract manifold M. are considered instead of 
phonon modes (so, there, the point of view involves a unifying framework for 
models of condensed matter physics) , such order parameters are assumed to be 
continuous across S. Here, on the contrary, we allow jumps of w. The proof 
below follows the one of the theorem in [dFM] quoted above. We adapt it to 
the situation envisaged here with slight modifications and report it with a cer- 
tain number of details just for the sake of completeness. In any case when in 
multifield theories the order parameter field takes values on a linear space, it 
can be considered discontinuous at the sharp discontinuity surface (if it exists) 
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and the counterpart of theorem above holds. On the contrary, when the man- 
ifold Ai of substructural morphologies does not coincide with a linear space, 
since A4 has finite dimension, it can be embedded isometrically in an appropri- 
ate linear space. However, the embedding itself becomes a prominent part of 
modeling. Although the isometric embedding is preferable because it preserves 
the quadratic part of the substructural kinetic energy (if it exists as in IIC), 
in fact, such an embedding is not unique (as non-isometric ones) and also not 
'rigid'. Then, the selection of the appropriate embedding (if necessary in the 
case of abstract order parameters) is not simple and general criteria suggested 
by physical instances seems to be not known. 

Proof. 

Conditions assuring the invariance of <p with respect to changes of observers 
and relabeling are given by j^cj) | Sl= o, S2 =o,s 3 =o = 0, with i = 1,2,3. They 
correspond to 

F T T • V s tr + N T S ■ V s tr+<9 m • (Vw) ro = 0, (56) 

T • V s v = 0, (57) 

3-6/ w «w))+§-V s £v w «w)) = 0. (58) 

If we shrink bs to bs (~l S uniformly in time, we get the pointwise balance 
(see [dFM] ) 

-[Q]U+ [5] -m + Dw s X = 0, (59) 

as a consequence of the arbitrariness of bs. 
If f 2 acts alone, then 

£ = T T v, Q = px-v, 3 = -P T v, (60) 

so that, as a consequence of (|57ll . Div^X = v-Div^- The arbitrariness of v 
and its continuity across E implies Ij51(l from 1)59(1 . 
If G acts alone, we get 

X = S T fv w «w» , Q = pw-£ Vw «w» , d = -S T 6/ W « w » » ( 61 ) 
and, from H58|l. 

DiwX = £ Vw «w» • (DivvS - i) ■ (62) 

Then, from l|59|l we obtain l|52fl thanks to the arbitrariness of the element £ 
selected in the Lie algebra of G. 
If f 1 acts alone, then 

Q = - i o F T x-ro-p(Vw) T w-ro, (63) 
(Qp |x| a + ip|w| a )l-P)» (64) 
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£ = -Cf an W-cv m , (65) 

with C tan and c defined respectively by ijSljl and lfS3|l and v m = tt> • m. 
Terms of equation l|59|) then become in this case 

-[Q]U + [31 • m = Po[F T xP ■ tu+p[(Vw) T w]U ■ to+ 



+ -p [|x| 2 ]ro • m + -p[|w|> • m - [P]tt> • m, (66) 

Div s (Cf an rt> + cv m ) = m ■ {Div s C tan + (Div^c) m) , (67) 

where the second equation is a consequence of Ij56(l . the circumstance that 
(I — m (g> m) • Vstt) = ((Vtr) m) • m, since tx> is isocoric, and properties 3 and 4 
of the definition of the relabeling f 1 of Bo including E. 

By inserting (|66|l and (|67|l in (|59|l . the arbitrariness of to implies 

p [F T x]C/ + p[(Vw) T w]C/ + ip [|x| 2 ]m+ 

+ i j o[|w| 2 ]m = [P T ]m + Dw s C tan + (D»«ec) m (68) 

and we shall evaluate the component along m of (|68|) . 

By indicating with v the averaged velocity v = (A) + U (F) m and using the 
relation [x] = — J7[F]m mentioned previously, we then get 

Po [F T ±]U ■ m = -~[p |x - v| 2 ] = Po [x] ■ v - p [ |x| 2 ]; (69) 

where \{p$ |x — v| 2 is the relative kinetic energy referred to S. We also get 

^ [|x| 2 ]^- Po [x].v+ip C/ 2 [|Fm| 2 ]. (70) 

by using once more [A] = — C/[F]m and the definition of v. By evaluating 
the normal component of (|68|l . using i|69[) . (|70() and taking into account that 
m'DwsCtan = Ctan ■ L, as it is simple to verify (see Lemma 2 in [dFM]), we get 
Ij53(l and the theorem is proven. 

Remark 6. Of course, for unstructured interfaces, i.e. in absence of surface 
energy, interfacial balances at items 1, 2 and 3 of Theorem 5 become respectively 

[P] m = -p [A] U, (71) 



m- [P]m =pU (Vw) w 



[S}m=-p[w]U, 
1 



iFml 



(72) 
(73) 
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6 Phason friction 



Non-conservative phenomena may occur in quasicrystals and involve different 
mechanisms such as (for example) viscous effects or plastic flows. Below we 
discuss just possible viscous effects due to 'internal' friction of phason nature. 
To fix ideas we restrict ourselves first to the case in which just friction of purely 
local nature exists. To account for it we may follow different ways. A classical 
one (see [BKMR] and [MR] for deep remarks about it) involves an integral 
Lagrange-d'Alambert principle. Preliminarily, we recall that balance equations 
lO and come from a variational principle of the form 

5 j L (j 1 (77) (X, t)) d(vol)Adt = Q. (74) 

B x[0,t] 

To account for possible internal friction of pure phason nature, we may then 
consider the following Lagrange-d'Alembert principle: 

5 I J C(j l (ri)(X,t)) d{vol) Adt j + / z v ■ 5w d {vol) A dt = 0, (75) 

\BoX[0,t] / B o X[0,t] 

with 

(a) z" = z* (F,w,Vw,w) G T^V U „ 

(b) z v • w > 0, Vw. 

The property (b) declares that the 'viscous' self-force z v is purely dissipative. 
In other words, one may say that formally (as it will be clarified further in next 
section) a balance of phason interactions of the form 

pw = z + DivS (76) 

still holds in non-conservative case, but with z, the self-force of phason nature, 
admitting an additive decomposition of the form z = z eq + z v , with z eq the part 
coming from thermodynamic equilibrium as in Corollary 2 and z v of purely 
dissipative nature. The abuse of notation between (|13J) and Ij76(l is rather neg- 
ligible because though the measures of interactions involved (namely phason 
stresses and self- forces) are placed within different thermodynamical settings, 
their nature is the same: they represent in fact interactions between neighboring 
material elements (the contact interactions represented by S) and interactions 
occurring within each material element. 

The inequality in (b) is satisfied by an expression of the type 

z" = cw, (77) 

with c = c (F, w,Vw, w) and c a scalar definite positive function such that 
c(F,w,Vw,0) = 0. In this case the balance of phason interactions becomes 

pw = Div (<9v w e) — <9 w e — cw. (78) 
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in the general case, while its reduced form for IQ is given by 

cw = Div (<9 Vw e) . (79) 

Of course, property (b) may imply explicit forms of z v more articulated than l|77|l 
and involving a structure of the type z v = Aw with A a second order definite 
positive tensor. Constitutive structures for z v involving tensor coefficients may 
also satisfy frame indifferent conditions (see discussions in [A] and [Si] about 
viscous stresses in classical viscoelasticity). 

Equation (|79|l fits the minimal model proposed in [RoLo]. 

We may consider also phason friction effects of weakly non-local (or better, 
gradient) nature. To this end we may consider not only a thermodynamic non 
equilibrium part z" of the internal self- force appearing in the sum z = z eq + z v 
but also a dissipative phason stress S v satisfying the decomposition S = S eq +S V 
with S eq as in Corollary 2 and S v of purely dissipative nature expressed by the 
inequality 

z v • w + S" ■ Vw > (80) 

that we presume to be satisfied by any choice of the rates involved. A possible 
solution of the inequality (|80|l is given by 

z v = c*w, S v = luVw, (81) 

with c* = c* (F, w,Vw, w,Vw), cj = u> (F, w,Vw, w,Vw), and c*, Co scalar defi- 
nite positive functions such that c* (F, w, Vw,0,0) = and u> (F, w, Vw,0,0) = 0. 
Of course, (|H 1|) is not the sole possible solution of (|80|l because tensor coeffi- 
cients may be involved and also linear combinations of w and Vw (see for a 
more general case [MA]). However, in the case of occurrence of (|80(1 . the reduced 
balance l|79|) for IQ becomes 

c*w = Div (<9 Vw e + wVw) . (82) 



7 SO (3) invariance and the nature of the bal- 
ance of phason interactions 

In deriving the balance (|13|l we have mixed the representation of phason in- 
teractions (obtained by means of the phason stress S and the phason self-force 
z) and their constitutive structure declared through the derivatives of the La- 
grangian with respect to Vw and w respectively. However, we have also adopted 
the formal counterpart of l|13fl . namely (|76fl . in non-conservative case, paying 
attention in the second circumstance to 'viscous'-like parts of the interactions. 
A question is whether such a pointwise balance holds formally always before 
discussing constitutive issues. Another connected question is whether an in- 
tegral (global) version of the pointwise balance of phason interactions can be 
postulated a-priori as a balance of phason-momentum. 

• The answer to the first question is affirmative: the balance of phason 
interactions holds in the form (|76|l independently of constitutive issues. 
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• As regards the second question, though an integral version of (|75|) can be 
in principle postulated because Y w is a linear space, it is not necessary 
because just the integral balance of standard forces and a non-standard 
balance of couples suffices to get pointwise balances. 

To prove previous statements we leave constitutive issues out of consideration 
and try to represent interactions just in their purely geometric form as objects 
power conjugated with the rates of the descriptors of the morphology of the 
body, namely phonon and phason degrees of freedom. 

Let b any part of Bq , i.e. any subset of Bq with non-null volume measure and 
the same regularity properties of Bq- We presume that the part in b interacts 
with the rest of the body and the external environment through interactions 
of bulk and contact nature, the latter exerted through the boundary db. The 
external power V£ (x, w) of all interactions over b, a linear functional over the 
space of rates x and w, is then given by 

Vl xt (x, w) = / (b • x + /3 • w) d (vol) + [ (Pn • x + Sn ■ w) d (area) . (83) 

Jb Jdb 

Here, b represents standard bulk forces and is decomposed as b = p Q h + h ln 
where b is the objective part which is coincident with the analogous b in Corol- 
lary 1 while b m is of pure inertial phonon nature. (3 is of pure incrtial phason 
nature (if phason inertia exists) while P and S are respectively the first Piola- 
Kirchhoff stress and the phason stress as in previous sections. Pn represents 
the traction developing power in the relative change of place of neighboring 
material elements at the boundary 9b imagining the phason activity frozen. 
<Sn pictures interactions developed across the boundary 9b between neighbor- 
ing material elements which do not change place but display different phason 
activity. As pointed out above, at each X we get P G Hom(T^Bo,T*B) and 
S eHom(T^B ,T^Y w ). 

We now require the invariance of (x, w) with respect to classical changes 
of observers ruled by SO (3). For such changes, the time parametrized family 
of automorphisms acting on the ambient space E 3 is the one of isometries so 
that, as usual, if x* is the value of the velocity x after the change of observer, 
we have 

x* =c(t)+q(t) A(x-x )+x, (84) 

where c (t) is the translational velocity, constant in space, Xo a point chosen ar- 
bitrarily and qA 6 so (3) at each t. Moreover, still for such changes of observers, 
SO (3) itself acts also over Y w and we indicate with w* the rate w measured 
after the change of observer, we get 

w* = w + q (i) A w. (85) 

Then, the requirement of invariance is 

7>r(x^*) = n ext (x,w) (86) 
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for any choice of translational c and rotational q velocities and for any part b 
(see [M] for a more general setting involving abstract morphological descriptors). 

The arbitrariness of c and q and their independence of space imply from 
H86|) the integral balances 

[bd(vol) + Pnd (area) = 0, (87) 

Jb Jdb 

/ ((x - x ) A b + w A 0) d (vol) + ((x - x ) A Pn + w A Sn) d ( area) = 0, 

Jb Jdb 



which are the standard integral balance of forces and a non-standard (due to the 
presence of the densities of phason interactions) integral balances of moments. 
They are the sole global conservation laws associated with the killing fields of 
the metric in the ambient space. 
The arbitrariness of b implies 

b + DivP = (89) 

from JgZJ and 

ePF T = w A (0+DivS) + (Vw) T S (90) 

from l|88l) . with e Ricci's alternating symbol. 

The inertial component of b, namely b m , and the explicit expression of (3 
can be identified by requiring that their power is the opposite of the rate of the 
kinetic energy, i.e. 

d f 

— {kinetic energy in b} - (b in ■ x + /3 ■ w) d (vol) = (91) 
dt J b 

for any choice of b and of the velocity fields. When sound-like modes appear in 
phason activity so that the kinetic energy in b is given by ~ J b (^p |x| 2 + p |w| 2 ^ d (v 
as in IIC, the arbitrariness of b and of the velocity fields implies 

b m = -p H, (3 = -pw, (92) 

In this way l|59"|l reduces formally to (JT^Jl. Moreover, from (|9~0"|) we get two 

information: 

1. at each X e Bo, the term ePF — (Vw) S is given by the cross product 
between w and an element of T^V W ~ R 3 that we indicate with z; 

2. z is just equal to DivS — pw so that we get l|76f) . 

Remark 7. In summary, when we represent interactions due to phason 
activity in quasiperiodic crystalline structures, an internal self-force arises a 
priori just as a consequence of requirements of SO (3) invariance of the power. 
Constitutive issues render explicit its structure as a function of state. In this 
way we find that the conservative part of z disappear for IQ because the relevant 
elastic energy does not depend on w while its dissipative part may play a role 
as in GHJ) (see [RoLo]). 
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8 Dissipative evolution of sharp interfaces 



We consider now the case in which the motion of X described by the vector 
field xb is not virtual, rather it is the real irreversible motion of X which can be 
e.g. identified with the boundary of an evolving defect. We assume that other 
irreversible phenomena like phason or gross friction do not occur. Dissipation 
is associated just with the evolution of X that we presume also to be coherent 
in the sense specified in Section 5. 

Wc then introduce a dissipative surface driving force fs along X. Coherence 
allows us to write the condition of dissipativity just on the normal component 
so that we require 

(fs • m) U < 0, (93) 

for any choice of U. We then modify (|50l) by adding a surface source term of 
the type 

fs • tvd (area) . (94) 

ti E ns 



Proposition 1. Let X be a structured surface with surface energy 
assume that during the dissipative evolution of X 



Let us 



d_ 

dt 



/ Qd (vol) + 3 ■ nd (area) 

J os J db-n 



X ■ nd (length) 



0(b s ns) 



fs • tod (area) = 



(95) 



b s ns 



for any part bs of Bo crossing X. // C and <j) are invariant in the sense of 
Theorem 5 and fi acts alone, the evolution of X along its normal is ruled by 

m [P] m + Ctan • L + Div^c+ 



-2*> U 



|Fm|- 



-pU 



(Vw) T w 



m p 



(96) 



where fy, is a positive driving coefficient such that /s = fy, (m, U). 



The proof follows Theorem 5 basically. In addition one may realize that the 
condition H93J1 implies that the normal component /s = fs • m be of the form 
/s = —fy,U with /s a positive coefficient. 

Remark 8. In absence of phason activity, (|96|l reduces to the balance 
equation describing the dissipative evolution of surfaces in simple bodies as 
discussed in [G] in presence of bulk deformation and surface energy. When bulk 
deformation is also absent, we may get the anisotropic motion by curvature. 
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